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Abstract. - 

Glasses have an excess number of low-frequency vibrational modes in comparison with 
most crystalline solids. We show that such a feature necessarily occurs in solids with low 
coordination. In particular, we analyze the density D{lo) of normal-mode frequencies lo and 
the nature of the low- frequency normal modes of a recently simulated system 0, comprised of 
weakly compressed spheres at zero temperature. We account for the observed a) convergence 
of D{ijj) toward a non-zero constant as the frequency goes to zero, b) appearance of a low- 
frequency cutoff LO* , and c) power- law increase of lu* with compression. We introduce a length 
scale r which characterizes the vibrational modes that appear at uo* . 



There is something universal and mystifying about the low-energy behavior of amorphous 
solids|21 El- In comparison to most crystals, amorphous sohds have a large excess number of 
low- frequency vibrational modes. In glasses these excitations are seen in the low-temperature 
specific heat as well as in the spectroscopy of the vibration modes in the terahertz range. These 
excitations affect the heat transport [3] and might well play an important role at the liquid- 
glass transition 0] . Nevertheless, little is understood about the cause of such excitations. 
Whereas in a crystal the vibrations are simply plane waves, in a glass, even at low angular 
frequency u>, they are much more complicated. In this Letter, we show that an excess density 
of vibrational states is a necessary feature of weakly-connected amorphous solids such as 
systems with repulsive, short-range interactions. Our analysis elucidates the cause and the 
peculiar nature of these low-frequency excitations. This gives a new approach for studying 
some of the ubiquitous phenomena found in glasses. 

A dramatic illustration of excess low-frequency vibrations was found in recent computer 
simulations |S1 of soft-spheres with repulsive, finite range potentials at zero temperature 
and zero applied shear stress. These simulations were carried out as a function of the packing 
fraction, (j) above the jamming threshold, (/)c, where the liquid acquires rigidity and becomes 
an amorphous solid O'Hern, Silbert et al. found that the average number of contacting 
neighbors per particle, z, the pressure, and the shear modulus vary as a power of [4> ~ 4>c)- 
Moreover, these simulations reveal unexpected features in the density of vibrational mode 
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Fig. 1 - D{uj) vs. angular frequency to for the simulation of Ref 1024 spheres interacting with 
repulsive harmonic potentials were compressed in a periodic cubic box to volume fraction (fi, slightly 
above the jamming threshold (/>c. Then the energy for arbitrary small displacements was calculated 
and the dynamical matrix inferred. The curve labeled a is at a relative volume fraction (f) — <j)c = 0.1. 
Proceeding to the left the curves have relative volume fractions 10~^, 10"'^, 10"'*, 10~*, respectively. 
Inset: Scaling of u)* vs. 5z. uj* for each {(p — 4>c) is determined from the data in the main panel as 
the frequency where D{uj) is half of the plateau value. 5z vs. (0 — (f>c) is obtained from the scaling 
measured in 0. The line has slope 1. 

frequencies, D{lu): (a) As shown in Figure 1, at </) = (jjc, when the system is most fragile, 
Dibj) has a plateau extending down to zero frequency with no sign of the Debye density 
of states normally expected for a three-dimensional solid, (b) As shown in the inset to that 
figure, the plateau erodes progressively below a frequency oj* that increases as 0— 0c increases, 
(c) The value of D{lo) in the plateau is unaffected by this compression. Similar behavior was 
recently seen in models of tetrahedrally coordinated covalent glasses [7]. Earlier simulations 
of a Lennard-Jones glass had also indicated an increase in D{uj) at low uj when z was lowered 
0. 

At the jamming transition, 0c, if particles with no contacts are excluded, the system 
network of contacts was found to be isostatic. That is, there are precisely the right number of 
inter-particle contacts to prevent free collective motions and the only zero-frequency modes 
are the uniform translations and rotations of the system as a whole. We shall see that if 
any contact were to be removed, one normal mode would become soft, with lu = 0. Soft 
modes have been discussed in relation to various weakly-connected networks such as covalent 
glasses iniEIl, Alexander's models of soft solids m, and models of static forces in granular 
packs[ni[TS|. 

In this paper we explain the plateau of excess low-frequency modes seen in Fig 1 by 
relating these modes to a certain set of soft modes. Our procedure for identifying the lowest 
frequency modes resembles that used for an ordinary solid. An isolated block of solid has 
three soft modes that are simply translations along the three co-ordinate axes. If the block 
is enclosed in a rigid container, translation is no longer a soft mode. However, one may find 
the lowest-frequency, fundamental modes by making a smooth, sinusoidal distortion of the 
original soft modes. We follow an analogous procedure to find the fundamental modes of our 
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isostatic system. First wc identify the soft modes associated with the boundary constraints 
by removing these constraints. Next wc find a smooth, sinusoidal distortion of these modes 
that aUows us to restore these constraints. 

The number of soft modes found in this way is far more than the three fundamental modes 
of an ordinary solid. For a d-dimensional cube of side L it is of order L'^"^ . As in an ordinary 
solid, the fundamental frequency uj^ of the distorted modes is of order L^^. The density of 
these modes D{ujl) (per unit volume L'^) is thus of order (L'^~^/L~^)L^^ ^ L°. In contrast 
to an ordinary solid, the density remains fixed even as L ^ oo and uj^ — > 0. 

In the following we make our procedure explicit for the simulation of Ref. PP. Using 
variational bounds, we show that D{lul) cannot go to zero at low frequencies. We extend our 
argument to set a lower bound on D{uj) for an extended range of frequencies. We account for 
the effect of compression in suppressing D{uj) at small lj. Finally we discuss the applicability 
of our mechanism to real glasses. 

Following we consider N soft spheres packed into a periodically-continued cube of 
side L at volume fraction (p. For inter-particle distance r < a, the particles are in contact 
and interact with a potential V{r) = |(1 — where a is the particle diameter and e a 
characteristic energy. For r > a the potential vanishes and particles do not interact. Particles 
with r < a have non-zero mutual energy and are said to be in contact. Henceforth we express 
distance, energy and mass in units of cr, e, and m, the particle mass, respectively. The energy 
SE relative to the equilibrium energy can be written as a quadratic form in the particle 
displacements from equilibrimir SRi...5Rn ■ 

+ J E {(^^^ - SR,).nrjy + 0{6R'} (1) 

where the sum is over all contacts (ij), r^J is the equilibrium distance between particles 
i and j, fiij is the unit vector along the direction ij, and {5Rj — 6Ri)^ indicates the projection 
of 5Rj — 6Ri on the plane orthogonal to n^. The term in brackets is proportional to the 
applied stress, and is much discussed in In continuous elastic media it is for example 

responsible for the fast vibrations of strings and drumheads. Near the jamming transition 
the contact forces vanish, r^J 1, and this term becomes arbitrarily small. In this letter we 
neglect it. 

It is convenient to express Eq ^ in matrix form, by defining the set of displacements 
6Ri...6Rn as a 3A^-component vector |(5R). Then Eq.(^) can be written in the form SE = 
(5R|A^|(5R). The corresponding matrix A4 is known as the dynamical matrix[Tl]. The 37V 
eigenvectors of the dynamical matrix are the normal modes of the particle system, and its 
eigenvalues are the squared angular frequencies of these modes. 

If the system does not have enough contacts, Ai has a set of modes of vanishing restoring 
force and thus vanishing vibrational frequency. For these soft modes the energy (the non- 
bracketed term of Eq.QJ in our approximation) must vanish. Thus the soft modes must 
satisfy the Nc equations: 

{SRi — 5Rj).nij ~ for all contacts (ij) (2) 

These linear equations define the vector space of displacement fields that conserve the distances 
to first order between particles in contact. Eq.® is purely geometrical and does not depend 
on the interaction potential. Each equation restricts the 3iV-dimensional space of |(SR) by one 
dimension. In general, these dimensions are independent, so that the number of independent 



SE = 



-1) 



[{SRj - SR,y 



Zr ■ 




Fig. 2 - Illustration of the boundary contact removal process described in the text. Eighteen particles 
are confined in a square box of side L periodically continued horizontally and vertically. An isostatic 
packing requires 33 contacts in this two-dimensional system. An arbitrarily drawn vertical line divides 
the system. A contact is removed wherever the line separates the contact from the center of a particle. 
Twenty-eight small triangles mark the intact contacts; removed contacts are shown by the five white 
circles. 

soft modes is 3A^ — Nc- Of these, six modes are dictated by the translational and rotational 
invariance of the energy function (5i?[(5R]. Apart from these, there are ZN — Nc — Q independent 
internal soft modes. 

A system of N repulsive particles is stable if it has no internal soft modes. This leads 
to the Maxwell criterion for rigidity A^^ > iN — 6. This is a non-local criterion for 
stability against collective motion of particles. At the jamming threshold, (pc, this inequality 
becomes an equality |S1 1121 113L 116) and the system is isostatic with coordination number 
Zc = 2Nc/N 6, as was found in [T] (i). 

In order to explain the large excess of low frequency modes in isostatic systems, we identify 
a large set of modes with low frequencies using a variational argument. Specifically, we shall 
perturb the boundary conditions in order to reveal a number p' ^ modes whose frequencies 
arc of order 1 / L. In order to show the existancc of this many modes, we construct a comparable 
number of variational trial modes. Because is a positive symmetric matrix, if an arbitrary 
normalized mode |(5R*) has 5E ~ {5Y{.*\M.\5Y{.*) = ujt^ we know that the lowest-frequency 
eigenmode has a frequency uji^ < ujt. Such an argument can be extended to a set of modes: if 
there are m orthonormal trial modes with ((5R*|A^|(5R*) < wt^, then there are at least m/2 
eigenmodes with frequency less than ^/2ujt (^) • 

For concreteness we consider the three-dimensional A^-particle system S of Ref with 
periodic boundary conditions at the jamming threshold. We label the axes of the cube by x, 
y, z. S is isostatic, so that the removal of n contacts allows exactly n displacement modes 

(^)This differs from rigidity percolation where contacts are deposited randomly on a lattice 1171 . In this case 
when a rigid cluster percolates, it contains over-constrained regions and is not isostatic. 

(^)If nia is the a'th lowest eigenvalue of M and if is an orthonormal basis such that {ea\M\ea) = 
then the variational bound of A. Horn [Am. J. Math 76 620 (1954)] shows that < ^^n^. Since 

V^v ^ X^i ^'I'i since ttIq > — (p/2)'"p/2i we have pUp > {p/2)mp/2 as claimed. 
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with no restoring force. Consider for example the system S' built from S by removing the 
p ^ contacts crossing an arbitrary plane orthogonal to (ox) (by convention at x = 0). 5', 
which has a free boundary condition instead of periodic ones along (ox), contains a space of 
soft modes of dimension p. (Fig. 2 illustrates how contacts are removed by changing periodic 
boundary conditions to free ones in a two-dimensional system.) In general when a contact (ij) 
is cut in an isostatic system, the corresponding soft mode is not localized near (ij), but extends 
over the entire system. This extended character comes from the non-locality of the isostatic 
condition that gives rise to the soft modes. This was confirmed in isostatic simulations jl2j. 
which observed that the amplitudes of the soft modes were spread over a non-zero fraction of 
the particles. Below we assume that a vector space encompassing a non-zero fraction of the 
soft modes have this extended character. We then consider the possibility of non-extended 
modes. 

We now use the vector space of dimension p' ^ of extended soft modes of S' to build 
p' orthonormal trial modes of S of frequency of the order 1/L. Let us define |(5R/3) to be a 
normalized basis of this space, 1 < /3 < p'. These modes are not soft in the jammed system S 
since they deform the previous p contacts located near a; = 0. Nevertheless a set of trial modes, 
|(5Rp, can still be formed by altering the soft modes so that they do not have an appreciable 
amplitude at the boundary where the contacts were severed. We seek to alter the soft mode 
to minimize the distortion at the severed contacts while minimizing the distortion elsewhere. 
Accordingly, for each soft mode f3 we define the corresponding trial-mode displacement (i|(5R*) 
to be: 

mSIi;)^C,sm{^){m,) (3) 

where the constants Cp are introduced to normalize the modes. The extended character of the 
soft modes implies that the Cp are of order one, and therefore bounded above by a constant c 
independent of L and /3 (^). The sine factor suppresses the gaps and overlaps at the p contacts 
near x = and x = L. Formally, the modulation by a sine is a linear mapping. This mapping 
is invertible if it is restricted to the extended soft modes. Consequently the basis |(5R/3) can 
always be chosen such that the |i5Rp are orthogonal. This procedure generates a 5E for each 
trial mode of order Indeed using Eq.QJ and Eq.(j2J), one obtains: 



<C^Tr/Lfz,nax=Ojl (5) 

where is the unit vector along (ox) , and Zmax is the maximum coordination number for one 
particle, z^ax = 12 for spheres (■*). Thus, we have accomplished our goal of finding ~ trial 
orthonormal modes below a frequency lul ^ We can now apply the above mentioned 

variational argument to show that the average density of states is bounded by a constant 
below frequencies of order ujl. 

(^)A normalized mode /3 is said to be extended if there is a constant C2 independent of N such that ((i|(5R^))'^ > 
C2/N for a non- vanishing fraction of the particles i. It is then straightforward to show that Cp defined as 
= Xy<ij> ^^'^'^(tt^^-'I'^'^/S)'^ bounded below by a positive constant ci independent of N. 

(*)In a polydisperse system Zmax could a priori be larger. Nevertheless Eq.Jl) is a sum on every contact where 
the displacement of only one of the two particles appears in each term of the sum. The corresponding particle 
can be chosen arbitrarily. It is convenient to choose the smallest particle of each contact. Thus when this sum 
on every contact is written as a sum on every particle to obtain Eq.JSJi the constant Zmax still corresponds to 
the monodisperse case, as a particle cannot have more contacts with particles larger than itself. 



6 



EUROPHYSICS LETTERS 



One may ask if the present variational argument can be improved, for example by con- 
sidering geometries of broken contacts different from the surface we considered up to now. 
When contacts are cut to create a vector space of extended soft modes, the soft modes must 
be modulated with a function that vanishes where the contacts are broken in order to obtain 
trial modes of low energy. On the one hand, cutting many contacts increases the number of 
trial modes. On the other hand, if too many contacts are broken, the modulating function 
must have many "nodes" where it vanishes. Consequently this function displays larger gra- 
dients and the energies of the trial modes increase. Cutting a surface appears to be the best 
compromise between these two opposing effects. Thus our argument gives a natural limit to 
the number of low-frequency states to be expected. 

We may extend this argument to show that the bound on the average density of states 
extends to a non-zero fraction of the modes of the system. If the cubic simulation box were 
divided into sub-cubes of size L /m, each sub-cube must have a density of states equal to 
the same {D{uj)) as was derived above, but extending to frequencies of order mujL. These 
subsystem modes must be present in the full system as well, therefore the bound on D{lo) 
extends to [0, ?tiljl]. We thus prove that the same bound on the average density of states 
holds down to sizes of the order of a few particles, corresponding to frequencies independent 
of L. We note that in d dimensions this argument may be repeated to yield a total number of 
modes, L'^~^, below a frequency lul ~ 1/-^, thus yielding a limiting non-zero density of states 
in any dimension. 

We now consider how these states are altered when the system is compressed to volume 
fractions above The simulations show that the extra-coordination number Sz = {z — z^) oc 
{(p — (/>c)°'^Em- Compression causes /S.Nc = N5z/2 ^ L^5z extra constraints to appear in 
Eq.(|21l. Cutting the boundaries of the system, as we did above, relaxes p ^ constraints. 
For a large system L^Sz > L^, thus p < ANc and Eq.ljSJ is still over-constrained so that no 
soft modes appear in the system. However, as the systems become smaller, the excess of ANc 
diminishes, and for L smaller than some /* ^ {6z)~^, p becomes larger than ANc and the 
system is again under-constrained. This allows one to build the trial modes described above 
in subsystems smaller than I*. These modes appear above a cut-off frequency uj* ~ (/*)^^; 
they are modes that contribute to the plateau in D(lo) above lo*. In other words, modes with 
characteristic length smaller than I* are little affected by the extra contacts, and the density 
of states is unperturbed above a frequency uj* ~ Sz ^^l- Using the explict potential V{r) 
defined above Eq.QJ ui* ~ (e/m)^/^cr~^5z. This scaling is checked numerically in the inset of 
Fig.l. It is in very good agreement with our prediction up to w 2. 

These arguments account for the occurrence of an excess in the low-frequency spectrum 
of D{u;) in amorphous systems of repulsive particles with finite-range interactions. The cor- 
responding low-frequency modes are related to the soft modes of free-boundary subsystems 
of size smaller than I*. I* does not appear in the static structure of these amorphous solids, 
but only in their responses. /* depends on the connectivity z and diverges at the jamming 
transition. At length scales larger than /*, we expect the system to behave as a continuous 
body. 

The argument above presumed that one could find a substantial subspace of p' extended 
modes in the space oip soft modes created by breaking the boundary contacts. This presump- 
tion, while plausible, is not trivial. Indeed, one can readily invent ways of breaking contacts 
that lead to most of the soft modes being localized near the broken contacts rather than being 
extended. For example if one breaks all the contacts in a spherical region, most of the spheres 
become unconstrained, and most of the modes are just the local motions of individual spheres. 
More generally compact arrangements of broken contacts lead to localized soft modes. We 
have checked numerically that the soft modes defined above Eq.Q are extended in the sense 
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needed 

The approximate low-frequency modes defined by Eq. ^ may help explain several puzzling 
properties of glasses. For example, their transport properties arc affected at temperatures that 
correspond to the frequency where the excess of vibrations appears |3| ■ Our analysis can give 
deeper insight into the vibrational modes responsible for e.g.. thermal conductivity. Moreover, 
our picture of these modes suggests a novel mechanism for mechanical and thermal relaxation. 
They give an alternative to the conventional cage picture in glasses: the non-locality of the 
soft modes suggests relaxation via extended events, such as the nonlinear buckling of modes 
resembling those of Eq.Q, rather than by local jumps out of a cage. The non- locality is a 
fundamental consequence of the Maxwell criterion for rigidity that underlies our argument. 

Finally we discuss the generality of this description of low-frequency vibrations of amor- 
phous solids. On the one hand, because our argument depends only on the connectivity of the 
system, it should be equally valid for systems with three-body forces (e.g. the model of Si02 
of [7| which leads to density of states similar to Fig.l) or particles with friction where degrees 
of freedom of rotation matter, provided the co-ordination number z is appropriately small. 
On the other hand, our approach neglects forces between non-adjacent particles, such as the 
Van der Waals forces present in all molecular glasses. Since these non-adjacent contacts are 
weak, our approach may be instructive in these real systems as well. 

We thank Leo Silbert for providing the numerical data and for discussions regarding the 
simulation results, and J.P Bouchaud, Bulbul Chakraborty, L.E. Chayes, P.G. de Gennes, 
Andrea Liu and Corey O'Hern for helpful discussions. We also acknowledge the support of 
the CFR fellowship, MRSEC DMR-0213745 and DOE grant DE-FG02-03ER46088. 
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